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The dynamical properties of the S = 1/2 antiferromagnetic XXZ chain are studied by the 
exact diagonalization and the recursion method of finite systems up to 24 sites. Two types 
of the exchange interaction are considered: one is the nearest-neighbor type, and the other 
is the inverse-square one. As the Ising anisotropy becomes larger, there appears a noticeable 
difference in the transverse component S xx (q, oj) between the two types of the exchange. For the 
nearest- neighbor type, the peak frequency of S xx (q, oj) for each q approaches the center of the 
continuum spectrum. On the contrary, the peak frequency for the inverse-square type moves to 
the upper edge of the continuum, and separates from the continuum for the anisotropy larger 
than the threshold value. Whether the interaction between domain walls (solitons) is absent or 
repulsive in the Ising limit leads to this difference in the behavior of S xx (q, oj). In the longitudinal 
component S zz (q,oj), on the other hand, the feature of the dynamics is scarcely different between 
the two types. The energy gap and the static properties are also discussed. 
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§1. Introduction 

The physics concerning one-dimensional (ID) quan- 
tum systems has been a long-standing problem. Dis- 
covery of many quasi-lD quantum magnets has spurred 
on the studies of the systems. Among those, there are 
5=1/2 systems in which a finite energy gap opens in 
the magnetic excitation spectrum. These systems have 
various origins of an excitation gap. For example, in 
CuGeC>3, the first inorganic compound in which a spin- 
Peierls transition was observed,!!! a gap appears below 
the critical temperature because of the lattice dimeriza- 
tion and the frustration. The other case with a gap is the 
XXZ chain with Ising anisotropy. CsCoCk-js known to 
be a typical ID Ising-like antiferromagnet.ePOne of the 
experimental means toonjeasure the spin gap is the inelas- 
tic neutron scattering.Ba' The neutron scattering brings 
us not only the information_nf the spin gap but also the 
whole excitation spectrumBEP In this context, the dy- 
namical structure factor (DSF) is a physical quantity to 
investigate carefully, because it is closely related with the 
experimental data detected by the neutron scattering. 

Besides connection with experiment, the study of dy- 
namics is indispensable in order to understand the na- 
ture of elementary excitations in the system. To this end, 
it is desirable to find the exact expression of the DSF. 
Haldane and ZirnbauerLr achieved this great work in the 
Haldane-Shastry (HS) modelJa'Er which is an isotropic 
spin model with exchange interaction proportional to 
the inverse square of the distance. Their result has the 
biggest merit that the picture of the elementary excita- 
tion is made simple. Recently, an integral representa- 
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tion of the exact DSF has b£en_clarified for the nearest- 
neighbor Heisenberg model JillEP although the result is 
considerably intricate. It is known that the HS model 
and the Heisenberg model belong to the same universal- 
ity class, and that the DSFs for these models are similar 
to each other. 

How do the features of the dynamics change when the 
Ising anisotropy is introduced in these models? This 
problem on the XXZ model is linked with formation 
of the finite gap. For the nearest-neighbor exchange, 
there exist a perturbation theory from the Ising limitp > 
and calculation based on the spin wave continuum. td-* 
To our knowledge, however, there is no work to study 
how the dynamics changes as the Ising anisotropy is in- 
creased. Moreover, for the inverse-square XXZ model, 
systematic studies of the dynamical properties have been 
lacking. 

The purpose of this paper is to clarify the effects of 
the Ising anisotropy and the interaction range on the 
DSFs for the ID XXZ model. Two types of the ex- 
change are considered: the nearest-neighbor exchange 
and the inverse-square one. We calculate the DSFs up 
to 24 sites via. the exact diagonalization and the recur- 
sion method. EfP We elucidate the change of the dynamics 
due to the Ising anisotropy and the difference between 
the two types of the exchange. Consequently, it is found 
that for the system with the large Ising anisotropy, the 
dynamics is sensitive to the range of the exchange in- 
teraction. In addition, the ground-state properties, the 
energy gap and the static structure factor (SSF) are in- 
vestigated. The energy gap and the SSF can be extracted 
from the process of calculation for the DSF. As a result, 
the characteristic difference is found in the opening of 
the energy gap between the two types of the exchange. 
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The remainder of this paper is organized as follows. In 
the next section (§||) we define the Hamiltonian of the 
ID XXZ model and the state vector. We review the re- 
cursion method in §[| Then, we show the results for the 
inverse-square XXZ model, comparing with those for 
the nearest-neighbor model in §[|. In §||, we discuss the 
behavior of the DSF from the Ising limit, and formulate 
the energy levels for the inverse-square model by using 
the two-soliton approximation. Finally we give conclud- 
ing remarks in 

§2. Model 

We consider the S — 1/2 ID antiferromagnetic XXZ 
model 



H • / '.- ( S ? S i + S i S d + AS i S . 



(2.1) 



where N is the number of spins which is assumed to 
be even, and denotes the exchange interaction be- 
tween spins i and j. The parameter A (> 1) represents 
the Ising anisotropy. The periodic boundary condition is 
imposed on the system. In this paper, we take two types 
of the exchange interaction: 

(i) the nearest-neighbor (NN) type, i.e., = JSjj+i; 

(ii) the inverse-square (IS) type, i.e., Jy = JD(xi — 
Xj)~ 2 , where D(xi — Xj) = (N/n) sin[7r(a;i — Xj)/N]. 
The Hamiltonians in the cases (i) and (ii) are written as 
Hnn and 7iig , respectively. 

States in the Hilbert space for a spin system can be 
represented by the position of up and down spins. If M 
denotes the number of down spins, the z-component of 
the total spin is given by 



JV 



^tot — E 



N 



M. 



(2.2) 



This quantity can be a good quantum number of the 
system. A state vector in the subspace with S* ot fixed is 
represented by 



M 



w = $>({*» n^W- 

{x} i=l 



(2.3) 



where \F) denotes the fully polarized up-spin state, and 
S~. is the spin- lowering operator at site a;,-. The wave 
function %p({x}) is symmetric in the position {x} = 
{xi, • • • , Xm} of the down spins. In eq. (|2.3|), the number 
of the bases is at Cm- Because the ground state of the 
system belongs to the subspace with 5f ot = for A > 1 , 
we can take the lowest state obtained by diagonalizing 
within the subspace S* ot = 0. The information of the 
ground state is necessary to calculate the DSF at zero 
temperature. 

§3. Numerical Method 

The DSF S ,1 ^(q,u>) is described as the Fourier 
transform of the dynamical correlation function 
(0\S£(t)Sjt +R \0), i.e., 

1 N N_1 r°° rlt 

1=1 R=Q 



5^ M^( g ,o; n ) 



(3.1) 



M""(g,w n ) = |W|0>| 

W n = E n — Eq, 



S'' 



(3.2) 
(3.3) 

= N-W^e-Msj* {q is a 
n) denotes an eigenstate of 



where /i = x, y or z 
momentum transfer), and 
Ti. with energy E n (Eq being the ground-state energy). 
S fJ,fJ, (q,uj) is_extracted from the Green function G w (g, z) 
as follows itS 

S^fau) = -- limlmG^(q,u + ir] + Eo), (3.4) 

G^(q,z) = (0\S» q (z-H)- 1 SZ\0). (3.5) 
is written in the form of a 



The Green function G MM (g, z 
continued fraction: 



G^(q,z) 



a 



(3.6) 



ai 



where S^{q) = (Q\St q S$\0) means the SSF. The coef- 
ficients a n and b n are obtained from the Lanczos algo- 
rithm in which we set the normalized initial vector as 
|/o> = ££|0>/||S£|0>||, 

l/n+l) - b-l x [(H - a n )\f n ) - ftnl/n-i)], 



a-n = {fn\H\f n ), 

b n+ i = \\{H- a„)\f n 



b n \fn-l 



b = 0.(3.7) 



Note that is normalized by b n+ \. 

The nonzero components of the DSF are S xx (q,u>), 
S yv (q,u>) and S zz (q,cu). Because the Hamiltonian eq. 
(2.1) has the axial symmetry, the transverse components 
S xx (q,u) and S vv (q,uj) are related with (q,u) and 
S~ + (q,ui) as follows: 



s xx {q,u) = syy{ q ,u) 



1 



^S+-(q,u;) 



2 b 



(3.8) 

In substance, we have only to calculate two quantities: 
the longitudinal component S zz (q, u>) and the transverse 
one S + ~(q,u) [= 2S xx (q, u)]. 

Thus the procedure for calculation of the DSF is sum- 
marized as follows: 

(1) We calculate the ground-state energy Eq and eigen- 
vector |0) via the Lanczos method and the inverse iter- 
ation one. In this case, we can take an arbitrary vector 
as the initial vector, except a vector orthogonal to the 
ground-state eigenvector. 

(2) We make |/ ) = S v q \0) /\\S v q \Q)\\ (y = z or -), and 
evaluate the coefficients a n and b n by using eq. (3.7). 

(3) Once we obtain a n and 6„, we can construct 
S zz {q,Lu) or S xx (q,u) from eqs. (3.4) and @. 

Originally the number of the coefficients a n and b n 
should be equal to the one of poles at each momentum q. 
For our systems, however, the results in the low-energy 
region are not sensitive to the number of coefficients used 
for n > 30. We calculate with use of at least 30 coeffi- 
cients a n and b n . 

In the practical numerical calculation, one needs to 
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take 77 as a small but finite value. This means that the 
(5-function is approximated by the Lorentzian. Thus 77 
corresponds to the Lorentzian width. For the purpose of 
inspecting the distribution of poles in detail, one should 
choose the value of rj as close as possible to zero. We use 
77 = O(10 -5 J) on taking account of computing time. 

Finally, we check the numerical precision by using the 
sum rule 

S^{q,u) 



N- 



duo- 



= 1. 



(3.9) 



All the data in this paper reach the precision more than 
0.999 for every momentum q. However, we note that the 
values of poles in the high-energy region are not reliable 
in contrast to the excellent aoquraac with the low-lying 
states, as the previous papersllStlJ'Llj pointed out. 

In this way, we can obtain uj n and M ^1 ' i (q , , u> n ) for 
each q in finite size N . In order to extract the behav- 
ior of S^^q, u>) in the thermodynamic limit from finite- 
size data_ we use the Lanczos spectra decoding (LSD) 
method.EZP The procedure gives the DSF as follows: 

M^(q,cj n+1 )+M^(q,uj n ) 



S^{q,w n ) = \ 



W„ + LOn+l 



-, (3.10) 



(3-11) 



However, this method is not easy to apply to cases where 
there exists an isolated mode with 5-function contribu- 
tion, and/or different kinds of continua overlap (for in- 
stance, a two-spinon continuum and a four-spinon one). 
If it is possible to separate each contribution from the 
raw data, we can use the remaining data to apply the 
LSD method. 

§4. Results for the Inverse-Square Exchange 
and Comparison with the Nearest-Neighbor 
Exchange 

We focus on the following physical quantities: the 
ground-state energy, the energy gap, the DSF and the 
SSF. For the IS exchange, the completely isotropic prob- 
lem (A = 1) has been solved analytically. After recalling 
the analytic expressions in each subsection, we present 
our results numerically obtained for A > 1 and discuss 
how these deviate from the ones for A = 1 . For compar- 
ison, we show the results for the NN XXZ model. Then 
we can bring out the features of the IS XXZ model in 
full relief. 

4-1 Ground- state properties 

In the IS isotropic case, the following simple form (i.e., 
the Jastrow-typa function) yields the exact ground-state 
wave functionrrcP 



N/2 N/2 

Its eigenenergy Eq is derived as follows: 

— = -^fi + A 

NJ 24 I N 2 



(4.1) 



(4.2) 



7T" 

24 



go-tropic model. 



-0.41123- ■■. (4.3) 
the ground-state energy 



-0.44314- 



(4.4) 



For the 
is known asE 

^ = -ln2 + i 
NJ 4 

This value is a little smaller than that for the IS isotropic 
model given by eq. (4.3). 

In the following we investigate the anisotropic case. 
First, let us look at the ground-state energy for the IS 
model with A > 1. We extrapolate the values in the 
thermodynamic limit by using the data from N =6 to 24 
and the following polynomial: 



T(N) 



r(oo) + ci-^ 



( 4 - 5 ) 

we show the values 



where T = E /{NJ). In Fig. 1 
evaluated for A = 2, 5 and 10 and interpolate those. 

For the NN model, Walker expanded the result of eq. 
(4.4) into the anisotropic case (A > 1).EP His result is 
equivalent to the following expression, .ss/hich was later 
obtained by des Cloizeaux and GaudinO 



= - sinh * 



,n=l 



(1 - tanhn*) 



- cosh*, 
4 



Eq(A) 
NJ 

1 1 — I 

(4.6) 

where A = cosh* and < * < +00. The limit A — ► 1 
(i.e., * — > 0) of eq. (4.6) coincides with eq. (4.4). In 
Fig. 1, we show the A-dependence of the ground-state 
energy by the thin solid line. This figure indicates that 
the difference of the ground-state energy between the IS 
exchange and the NN one widens as the Ising anisotropy 
becomes larger. 

Next, we consider the wave function. Because eqs. 



(4.4) and (4.6) for the NN model are derived on the basis 
of the Bethe ansatz, the corresponding wave functions 
are very complex. From here, we concentrate the wave 
function for the IS XXZ model. 

In this case, the Jastrow-type function 



M 



M 



v J ({x})=n(-i) x< iii i5 ( a 



(4.7) 



i<3 



is the exact ground-state wave function in the subspace 
of the number M of down spins under the conditions 



A(M - 1) < N, 



(4.8) 



and A =even integer, where A = A(A — l)/ 2& When 
A > 1 (i.e., A > 2), this wave function eq. fl4.7|) is not 
the absolute ground state (i.e., M = N/2) of the system. 
In order to test the accuracy of this wave function as a 
trial function, the overlap |(V'exact|'0j)| between eq. (4.7) 
and the exact eigenfunction ?/>oxact is calculated by the 
numerical diagonalizaticp. 

In our previous paper ,£3 we showed that the overlap is 
larger than 0.95993 (of course, smaller than unity) over 
whole range of the Ising anisotropy (A > 1) for N = 12 
and M = 6. We extend this study up to N — 22. Figure 
2 shows the overlap against the Ising anisotropy for 22 
sites. In this case, the overlap is not as close to unity as 
that for 12 sites, but is still more than 0.88305. The A- 
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dependence of the overlap has a conspicuous minimum at 
a certain value of A. This value of A sweeps toward the 
smaller one as the size N becomes larger. For N = 6, 12 
and 22, the integer values of A with the minimum overlap 
are A =8, 6 and 5, respectively. 

One may naturally ask the size dependence of the over- 
lap for fixed Ising anisotropy. We show the cases A = 2, 5 
and 10 in Fig. 3. In all cases, the overlaps decay with 
increasing size N. However, the decaying rate is differ- 
ent: it is most rapid for A = 5. Fo r A = 1 and A — > oo, 
the Jastrow-type function eq. (4.7) represents the exact 
eigenstate, as mentioned in ref. 22. The results in Fig. 3 
suggest that the overlaps for 1 < A < oo vanish in the 
limit N — ► oo. Nevertheless the wave function eq. (4.7) 
is surely an excellent trial function for the finite systems, 
at least with TV < 22. 



4-2 Energy gap 

To compare with new results for the IS model, first we 
review the result for the NN model. Des Cloizeaux and 
Gaudin investigated the nature of the ground state and 
the first excited states for all values of Ao Particularly, 
they demonstrated that for A > 1, an energy gap Eq(A) 
appears, which can be defined as, 

E G (A) = lim [E(A, l,q) - E(A, 0, q)} . (4.9) 

Here E(A, S* ot ,q) means the lowest energy in the 
subspace with two quantities (S* ot ,q) fixed for the 
anisotropy parameter A. The expression of the energy 
gap is given by 



7rJsinh\l/ \ - 
E G (A) = £ 

n=— oo 



i 



cosh[(2n + l)7r 2 /(2*)]' 

(4.10) 

where A = cosh fy. This gap increases very slowly at the 
threshold as follows: 



Eg(A) ~ 47rJexp 



' v / A — I 



< A - 1 « 1. 



(4.11) 



Therefore the energy gap has the essential singularity at 
the critical point A = 1; all differentials become zero at 
the point. When A goes to infinity, the gap behaves as 



£ G (A)~ J(A-2). 



(4.12) 



In Fig. 4, we show En M.S M function of A by the thin 
solid line. 

Next, we numerically investigate the behavior of the 
gap for the IS XX Z model. The excitation gap is defined 

as, 



E q G (A)= lim [E(A, 1, q* 



)--E(A,0,0)]. (4.13) 



We consider the cases q* = and tt. 

In Figs. 5(a) and 5(b), we show the results of finite-size 
anal ysis with the data of N — 6-24 and the polynomial 
eq. (gj) [r = E(A, l,g*)-£'(A,0,0)]. In some cases the 
extrapolated value is slightly negative. For instance, the 
value of E G is -0.030023 for A = 2. These cases with a 
negative gap indicate the degree of numerical inaccuracy. 



The A-dependence of the gap is shown in Fig. 4. 

We find the following features from Fig. 4. First, for 
the IS type, the opening of the gap is slower than that for 
the NN one. Second, both E G and EZ for the IS model 



seem to be well fitted by the essentially singular form.fci$ 
However, it is difficult to evaluate the quantitative form 
by our numerical calculation: we cannot determine the 
accurate critical point. 

4-3 Excitation spectra and dynamical structure factors 
In the IS isotropic case (i.e., the HS model), Haldane 
and ZirnbauerEP derived the exact expression of the DSF 
in the thermodynamic limit. It is given by 



S xx {q,io) = S z 

ie(wu(«)- 

4 



w)0(w 



[{u-^{q)){u-^ + {q))YI* 
where 0(w) is the step function, and 



^L+(q) = -{q- tt){2tt - q), 

uxs(q) = ^i( 2lT - <?)• 



,(4.14) 

(4.15) 
(4.16) 
(4.17) 



Equations (4.15) and (4.16) represent the lower edge of 
the excitation spectrum for < q < tt and for tt < q < 2tt 
respectively, and eq. (4.17) gives its upper boundary for 
< q < 2ir. The intensity diverges at the lower edge of 
the continuum spectrum and jumps to zero at the upper 
edge. 

Here we perform the numerical calculation of the DSF 
in the HS roadel. The same calculation has been done up 
to 16 sites£-3> We show our result for N = 24 in Fig. 6. 
It reveals clearly the compact support for the excitation; 
there is no intensity outside of the continuum.E§> This 
is due to the selection rule of the Yangian, which is a 
special symmetry with the HS model. Namclx the DSF 
consists only of the two-spinon excitationsE'ES 

In this context, the situation is more complicated for 
the NN isotropic model (i.e., the Heisenberg model). For 
this model, the higher-order-spinon processes jas_wBll as 
the two-spinon ones contribute to the DSfJB'BEZP The 
exact two-spinon DSF has been obtained by an integral 
representation.^ j—. 

We test the LSD methodlilP in the HS model. In Fig. 
7, we compare the numerical result (closed symbols) pro- 
cessed by the LSD method for N =16-24 with the exact 
one (solid line) of eq. (4.14) at q — tt, which becomes 



S xx {ir,LU 

1/(4") 




for < uj/J < tt 2 /4, 
otherwise. 



(4.18) 



It is found from Fig. 7 that the numerical result agrees 
fairly well with the exact one, especially near the upper 
boundary. There appears clearly the step-function-likc 
behavior peculiar to the HS model. The numerical re- 
sult deviates somewhat from the exact one in the lower- 
energy region where the number of poles is fewer (that 
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is, the density of states is smaller). The deviation oc- 
curs because the finite diffcrenpc is substituted for the 
differential in the LSD methodO' 

Now, let us turn to the anisotropic case. We show 
S xx (q,u) for the IS XXZ model with A = 2,5, and 10 
in Figs. 8(a), 8(b) and 8(c), respectively. Characteristic 
features with increasing Ising anisotropy are: 

[A] The shape of the continuum tends to be symmetric 
about q/n = 1/2 but the intensity does not. 

[B] The main intensity shifts from the low-energy re- 
gion to the high-energy one of the continuum. 

Encouraged by the validity of the LSD method for the 
HS model, we apply the method to anisotropic cases. 
Unlike for A = 1, first we select poles with main in- 
tensity (> 10~ 2 ), which are presumably due to two- 
soliton contribution, and subsequently process the se- 
lected data by the LSD method. Figure 7 shows the 
results at q = tt for N =16-24 and A = 1, 2 and 5. The 
data for (N, A) = (18, 2) have been removed since some 
numerical inaccuracy causes the irregularity of the de- 
pendence on N. This figure makes the above feature [B] 
clear. Note that the peak intensity does not shift gradu- 
ally inside the continuum; it appears either at the lower 
boundary or at the upper one. 

For comparison, we apply the same method to the NN 
model. At q — tt, the number of poles with main in- 
tensity we have selected is equal to the number of all 
poles for the HS model in each N. In Fig. 9 we show 
the change of S xx (ir,u>) as a function of A. For A = 1, 
S xx Jjx, u>) has a rounded shoulder at the upper bound- 
aryEZI' in contrast to the HS model. Moreover, unlike 
the IS case, S xx (tt,uj) which forms the continuum be- 
comes roundish with increasing A. The peak frequency 
approaches the center of the continuum. 

What happens in the behavior of the DSF when the 
Ising anisotropy is much larger in both the IS XXZ 
model and the NN one? We show the result for the 
IS model with A = 10 in Fig. 10. To see the size depen- 
dence of poles and intensities, we show the behavior of 
M xx (tt, u>) [see eq. (3.2)] in the main plot of Fig. 10(a). 
The following features are observed: 

• Among poles with main intensity (which all are in- 
dicated in this figure for A = 10), the highest pole 
position is hardly dependent on size N, in contrast 
with other poles. 

• The intensity of the highest pole oscillates with N 
but seems to converge to a certain value. 

From these features, it is likely that the highest pole 
is an isolated mode separated from the continuum to 
which the other poles belong. For the purpose of judg- 
ing the separation, it is useful to see the behavior of 
NM xx (tt, oj), shown in the inset of Fig. 10(a). This fig- 
ure strengthens the possibility of the isolated pole. We 
regard the highest pole as the <5-function contribution. 
Then we eliminate this pole and analyze remaining poles 
by the LSD method, which results in the continuum part 
of S xx (n,u>) shown in Fig. 10(b). It seems to be a mono- 
tonically increasing function of to. On the other hand, 
there appears no isolated mode for the NN model with 
A = 10, which is consistent with the result by the per- 
turbation approach from the Ising limit x3 



Because there exists the isolated mode for the IS model 
with A = 10, one may naturally ask the threshold value 
A c at which the mode begins to appear. We judge that 
the isolated mode exists if the behavior of M xx (tt, ui) [or 
NM xx (ir, ui)] satisfies the following feature in addition 
to the two itemized above: two-soliton contribution with 
main intensity does not overlap in higher-order-soliton 
one with very small intensity. See Figs. 11(a) and 11(b), 
where we show the results for A = 6 and 7 respectively. 
We infer 6 < A c ^ 7 from these figures. In order to deter- 
mine the precise value of A c , it is necessary to perform 
detailed analyses of poles and intensities in systems with 
still larger size. 

Next, we show S zz (q,u) for the IS model with A = 
2 and 5 in Figs. 12(a) and 12(b), respectively. For 
S zz {q,u>), the character of the dynamics is scarcely dif- 
ferent from that for the NN model. Instead of [B], the 
peculiarity to S zz (q,ui) comes out: [B'] The intensity at 
oj ~ and q — tt is conspicuously enhanced as the Ising 
anisotropy becomes larger. It is because the Neel-type 
correlation strengthens at q — tt. Meanwhile, the other 
intensities are suppressed. For A = 5, the poles with 
weak intensity seem to belong to a continuum. This con- 
tinuum may have the feature [A] on the. analogy of the 
result obtained by Ishimura and Shiba.liif 

4-4 Static structure factors 

Before the IS isotropic model (i.e_ the HS model) was 
proposed, Gebhard and VollhardtEEf calculated the ex- 
pression of the SSF in terms of the Gutzwiller wave func- 
tionc$ as follows: 



S xx {q) = S* 



'( 9 ) = -iln(l-£ 

4 V TT 



(4.19) 



for < q < tt. After that, it was demonstrated that the 
ground-state wave function of the HS model is identical 
to the Gutzwiller one.S 0> Thus eq. (JOJ|) turns out to 
be the exact expression of the SSF for the HS model. 
This was confirmed by carrying out the integration of eq. 
(4.14) over W.EP From eq. ( 4.19| ), the SSF is proportional 
to q for small q, and it is logarithmically divergent at 
q = tt. 

In Fig. 13, we compare this exact result (solid line) 
with our numerical data (closed symbols) for the HS 
model (A = 1). The fit of the symbols to eq. Q4.19| ) 
is very good except at q = tt. The data at q — tt are 
noticeably dependent on the size N, reflecting the loga- 
rithmic divergence. 

Let us turn now to the anisotropic case. First, we look 
at the transverse component S xx (q). Figure 13 shows 
S xx (q) for the IS model with A = 2 and 5 besides A = 
1. The ^-dependence of S xx (q) becomes weaker with 
increasing Ising anisotropy. S xx (q) is independent of q 
in the Ising limit; the value of S xx (q) is equal to 1/4. 

Next, we show the longitudinal component S zz (q) for 
the IS model with A = 1,2 and 5 in Fig. 14. As the 
Ising anisotropy becomes larger, S zz (q) at q — tt rapidly 
increases; S zz (tt) has considerable dependence on size N. 
On the other hand, S zz (q) for q < tt are suppressed. This 
leads to the weak intensity in S zz (q,Lo) for < q < tt. 
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§5. Interpretation of S xx (q,ui) from the Ising 
Limit 

5.1 Qualitative interpretation 

What is the origin of different behavior of S xx (q, u) 
for A > 1 between the IS XX Z model and the NN one? 
We consider this question from the Ising limit for both 
models. 

First, let us take notice of the energy levels. In the 
Ising limit, the ground-state level is doubly degenerate 
corresponding to the Neel states with a different choice 
of sublattice for both models (see ref. 30 and Appendix) . 
However, the nature of the next level differs between 
the two models. Here we limit the Hilbert space to the 
subspace with S* ot — ±1, in order to inspect the levels 
related with S xx (q,oj). For the NN case, the next level is 
separated from Neel states by J and is highly degenerate. 
The states with this level have only two domain walls, or 
solitons. We can interpret as free solitons. On the other 
hand, for the IS type, the degeneracy is removed except 
that due to the translational symmetry. We find numer- 
ically that the first excited state is the state which has 
two solitons with the longest distance; moreover, among 
the excited states with two solitons, the state with the 
nearest-neighbor solitons has the highest energy (see Ap- 
pendix). This means that two solitons are repulsive for 
the IS Ising model. 

Next, we consider the transverse component S xx (q, u>) 
of the DSF. See the results of Appendix. In both models, 
the states contributing to S xx (q,ui) are limited to those 
with the nearest-neighbor solitons. As mentioned above, 
however, the energy level of this state is equal to that of 
other states with two solitons for the NN Ising model, 
whereas it is the highest among the energy levels of the 
two-soliton states for the IS model. 

When the spin-flip terms are introduced in the Hamil- 
tonian, the energy level is split up and down for the NN 
model. On the contrary the levels for the IS model are 
only modified a little because these are not degenerate 
from the beginning. The DSFs come to contribute to 
many states in both models. Nevertheless each nature 
in the Ising limit remains if the spin-flip terms are much 
smaller than the Ising term: most of the intensity comes 
from the vicinity of the energy level with the nearest- 
neighbor-soliton state in the Ising limit. 

In this way, the behavior of S xx (q, ui) for A ^> 1 is 
qualitatively explained. 

5.2 Two-soliton problem in the inverse- square XX Z 
model 

In order to treat the discussion on the energy levels 
quantitatively, we study the excitation spectrum for the 
IS XXZ model by using the two-soliton approximation. 
This approximation should be valid for description of the 
low- lying excitation near the Ising limit. 

We redefine the Hamiltonian as follows: 



Wis/A = Wig 



N 

E J > 

i<j 



i J 



Sa S ' ■ 



where J;. 



JD{xi 



(5.1) 

The two-soliton eigenstate 



is written by a linear combination of all possible con- 
figurations of two domain walls from the reference state 
|Neeli) = |UU"-U),i.e. J 

H>) =X>S^|Neeli> 

i 

+ E - ^)/*A + Sr+i ■ • ■ S+lNeeb) 

+ e( Xl - x,)f l3 S+Sr +1 ■ ■ ■ S+Sf ■ ■ ■ S+lNeeU)] , 

(5.2) 

where Q(x) is the step function. Note that both sub- 
scripts i and j of /y run over only even numbers. Thus 
the total degree of freedom becomes (N/2) 2 . 
Let us consider the eigenvalue problem 



H is \i>) = E\i/>). 



(5.3) 



First, we define a relative coordinate r between the two 
solitons as follows: 



(xj + N) — Xi 



for Xj > Xi, 
for Xj < x^ 



(5.4) 



Then the permissible values of f are limited to 
0,2,4, N — 2. By taking the inner product 
of both sides of eq. (5__3) with every basis state 
(Neeli|5~ • • • S' i * fl S'~, we obtain the following equation 
for the amplitudes /y: 



[E-Eo- E, 



add(f+l)] fij 
+ (1 — ^i,i-2 — ^j.i+JV-a) (/i,i+2 

where 



i-2,j) 



En 



48 



-JN 1 



2 

A2 



E. 



add(r+l) 



2 \NJ 



r+l N ( 

-EE - ! 



i=l j=f+2 
2 



sin [n(xi — Xj)/N] 
1 



= 0,(5.5) 
(5.6) 

(5.7) 



NJ mC{ir/N) 



(5.8) 



Next, we separate out the center of mass motion. Be- 
side the definition of the relative coordinate f by eq. 



(5.4), we introduce the center of gravity R, the total 
momentum q and relative momentum fc, 

for Xj > Xi, gs 
for Xj < Xi, 

(5.10) 



1R " Xi + {xj + N) 
q = fci + k 2 , 2k — k 2 — k\ 



where k\ and k 2 denote the individual momenta of the 
two solitons. Then the amplitudes can be expressed 
as follows: 



A 



ikr 



(5-11) 

where each summation is over < q, k < tt. This is be- 
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cause, both R and f go by two steps on the lattice spac- 
ingElP Substituting eq. ( 5.11 ) into eq. (5.5), we obtain 
the equation for F q (f): 

[E - E - £ ad d(f+i)] F g (f) - eJ(l)cos 9 

x [(1 - 5f,o) F q (f - 2) + (1 - Sr, N - 2 ) F q (f + 2)] = 0. 

(5.12) 

Since we consider the anisotropy near the Ising limit 
(i.e., e 1), the excitation energy to of the Hamiltonian 
eq. ([O]) is almost E — E where E is given by eq. (5.6). 
Furthermore, we multiply eq. (5.12) by exp(— ikf) and 
sum on all f, to obtain the equation for f q {k): 

{cj-u> q (k;N)]f q (k) = ^ J2 U(k,k')f q (k'), (5.13) 



where 

uj q {k;N)=E c + 2eJ(l) ( 1 



N 



add(f+l)) 



cosgcos(2fc),(5.14) 



(5.15) 



C/(fc,fc')=^c- i ( fe - fc ')^ a dd (F+ i) 

f 

-eJ(l)cosg (e" 2ifc ' 



+ e- 



2i/,- 



(5.16) 



The scattering states (i.e., two-soliton contin uum) ace 
obtained by ignoring the right-hand side of eq. (5J_3).c3 
The continuum is given by 

Co q {k;oo) = J [0.693147 + 2e cos q cos(2fe)] . (5.17) 

Here we have used 

J(l) -» J, (5.18) 

E c 0.693147 J, (5.19) 

for N — » oo. (The value of E c is the extrapolated one 
from the numerical calculation for the finite size.) 

Let us sec the validity of the two-soliton approximation 
for N = 24 and e = 0.1 (i.e., A = 10). In Fig. 15(a) we 
show the r esult obtained by solving the eigenvalue equa- 
tion ( 5.13 ) numerically. This result is compared with the 
one of the pole positions via the recursion method in Fig. 
15(b). This figure shows the poles with main intensity 
from two-soliton contribution, and exhibits together the 
ones shifted from [7r,27r] to [0,7r]. The numbers of the 
energy levels are the same between the two figures at 
each momentum q, except at q — 5n/12, ir/2 and 7w/12. 
In the two-soliton approximation, all of the degree of 
freedom appear within the region < q < n, i.e., half 
the Brillouin zone. We judge that this approximation is 
rather good for e = 0.1. 

Figure 15(a) suggests the existence of the isolated 
mode above the continuum spectrum. In order to make 
this clear, we show the solutions of eq. (5.13) for TV = 100 
and e = 0.1 in Fig. 16. First, we verify the good fit by 



the expression eq. (5.17) for the continuum spectrum. 
Next, one mode separates obviously from the continuum 
over the whole range of the total momentum q. In ad- 



dition, another mode is seen for 0.45^g/7r^0.55. These 
isolated modes are interpreted as the sets of the anti- 
bound states due to the repulsive solitons. The upper of 
the two modes corresponds to the isolated mode for the 
IS XXZ model with A — 10 obtained by the recursion 
method. 

Here we compare our result with the one for-.the NN 
XXZ model by using the same approximation£$ In this 
case, the two-soliton continuum is given by 



u>q(k, oo) = J [1 + 2ecos gcos(2fc)] 



(5.20) 



with < g, k < 7r in our definition. No anti-bound state 
appears when the corresponding eigenvalue equation is 
solved numerically. The reason is that the solitons are 
identified with free particles except for the hard-core re- 



pulsion, as mentioned in §5.1 



The repulsion between solitons turns out to arise from 
the interaction Jy with j > i + 2 in the Ising term of 
eq. ( jj~ll ) . The necessary interaction for the repulsion be- 
tween nearest- neighbor solitons is the oae. which reaches 
at least the next-nearest-neighbor spinS Independent 
of the range of the interaction, the terms with e in the 
Hamiltonian play two roles: one is to transfer the soli- 
tons (i.e., kinetic energy), and the other is to change the 
number of solitons. However, the latter role is neglected 
within the two-soliton approximation. The kinetic en- 
ergy causes the dispersion in the excitation spectrum. 

§6. Conclusion 

We have investigated the ground-state properties, the 
energy gap, the DSFs, and the SSFs of the XXZ chain 
with Ising anisotropy for the two types of the exchange 
interaction. The main results of this study can be sum- 
marized as follows: 

(1) For the inverse-square type, the overlap between 
the exact wave function and the Jastrow-type one tends 
to zero in the thermodynamic limit over the whole range 
of the Ising anisotropy, except the isotropic case and the 
Ising limit. 

(2) In each type, the dependence of the energy gap 
on the anisotropy parameter has the essentially singular 
form at the critical point. For the inverse-square type, 
the opening of the gap is slower than that for the nearest- 
neighbor type. 

(3) With increasing Ising anisotropy, the peak position 
of S xx (q, w) for the nearest-neighbor type approaches the 
center of the continuum; on the other hand, the peak po- 
sition for the inverse-square type moves to the upper edge 
of the continuum, and separates from the continuum for 
the anisotropy larger than the threshold value. 

On the basis of the result (3), we have proposed the 
picture of the repulsive solitons for the inverse-square 
XXZ model when the Ising anisotropy is large. The 
anti-bound mode in the (q, w)-plane of S xx (q,ui) is due 
to the repulsive solitons. In a quasi- ID Ising- like anti- 
ferromagnet where the exchange reaches a long range, 
the repulsion will lead to a neutron scattering cross sec- 
tion which is strongly enhanced at the upper boundary 
or above the continuum. In this connection, we mention 
that there are both two-magnon continuum and bound 
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mode for the alternating spin chainH The existence of 
the bound mode contrasts with the anti-bound mode in 
our system. 
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Appendix: Dynamical Structure Factors of the 
Inverse- Square Ising Model 

The Hamiltonians His an d Hnn are given in the Ising 
limit as follows: 



Wis/A 



Wnn/A 



JV 



1 



i<3 



N 



A- 



</E] SI S i+1 . 



SI Si , (AT) 



(A-2) 



In this Appendix, we calculate the DSFs for the IS Ising 
model defined by eq. (AT). 

It is not trivial that the ground state for the IS Ising 
model is theJNeel state. To this end, we quote Hub- 
bard's work.EP He demonstrated that for the ID model 
only with the interaction of electrons (negligent of the 
hopping term), the electrons of the ground state have a 
periodic arrangement which is regarded as the general- 
ized Wigner lattice under certain conditions. By adapt- 
ing his conclusion to the spin system, the following is 
found: the exact ground state of an Ising model is the 
Neel state, provided that the interactions J(r) satisfy the 
two conditions: 

J(r) -> as r -> oo, (A-3) 

J(r + 1) + J(r- 1) > 2J(r) for all r > 1, (A-4) 

where r denotes the distance between two sites. The IS 
Ising model fulfills this conditions for N — ► oo. 

Next, we calculate the eigenenergies of the Neel state 
and states with two domain walls. We begin with the 
Neel state. Let us write the two degenerate Neel states 
as iNeeb) = | fltl ■ ■ ■ Tl) and |Neel 2 ) = | HIT • • ■ ID- 
Then 

N 



Hisi|Necl 1 ) = lj(^) 2 ^ 1 
2 \N / *— -f sm \Tr(xi - 



i)/N] 



xi(-l) Xi -^|Neeli). 



(A-5) 



Using the following relation: 

E • / /ft n =- 7 (^ + 2), (A-6) 
^-J sm (-nr/N) 6 

the eigenenergy Eq (i.e., the ground-state energy) is 



given by 



Eo = — 77: JN[ 1 + — 7T 



2 

N 2 



(A-7) 



The eigenenergy of iNeekj) is also equal to Eq. Further- 
more, we define the state with two domain walls as 

lDw(m)> = 1 n_^m u • • • ut) , (a-s) 

where m is the distance between two domain walls and 
should be odd in the subspace S* ot = 1 (or —1). Then the 
eigenenergy i£ D w(m) of the state |DW(m)) is obtained as 
follows: 



^DW(m) — E { 



1 



-J 



N 



E E 



~! J= m+i sin t 71- ^ _ x j)I n \ 

(A-9) 

Particularly for m = 1, i.e., the nearest-neighbor do- 
main walls, the double summation of eq. ( |A-9| ) can be 
performed to give £"dw(i) = Eq + Equw with 



2 



(ATO) 



Also by calculating eq. ( |A-9[ ) numerically, it is found that 
Bdw(ih) is a monotonically decreasing function of m in 
the range 1 < m < N/2. 

Now that we have found the necessary eigenenergies, 
it is easy to calculate the DSFs for the IS Ising model. If 
|GSi) = (|Neeli) - |Neel 2 ))/\/2 and |GS 2 ) = flNeeli) + 
|Neel2))/\/2, the longitudinal component S zz (q, to) is de- 
rived as 

2 



d=l n 



(ATI) 



This is the same as the result for the NN Ising model. 
On the other hand, the transverse component S xx (q,u) 
is calculated as follows: 



1 



"few) 



I E E I |GS d ) | 2 <5 (c - (S„ - So)) 

(AT 2) 



d=l n 



4 ^ 

d=l 



M d S (w - (E n - E )) , 



where 



M, 



-Y 



N 



Ee~ i9 ^7|Neel a 



(ATS) 



The states \n) with nonzero matrix element are limited 
to those with only one flipped spin from the Neel state, 
i.e., |DW(1)) and several of its translationally invariant 
states; the total number of these states is N/2, so that 
Mg = 1/2 (d = 1,2). In addition, all these states have 
the eigenenergy -E'dw(i) = Eq + Equw- As a result, the 
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following equation is obtained: 

1 



Riera and B. C. Sales: preprint (cond-mat/9708078) 



S xx {q,u) = -5(lo-E gbw ). 



(A-14) 



Incidentally, S xx (q, uj) = (1/4) -S(u- J) for the NN Ising 
model given by eq. (A- 2). 
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Fig. 1. Ground-state energy per site as a f unct ion of A. The thin 
solid line is the exact result given by eq. (4.6) for the NN XXZ 
model. The open circles show the data extrapolated from finite 
size data (jV = 6-24) for the IS XXZ model with the second- 
order polynomial fit. The value for A = 1 corresponds to the 
exact solution given by eq. (4.3) of the HS model. The thick line 
is guide to the eye. 



Fig. 2. The overlap vs the Ising anisotropy A for the IS XXZ 
model. N and M denote the numbers of sites and down spins, 
respectively. 



Fig. 3. Size dependence of the overlap for various values of A. 



Fig. 4. Energy gap as a functi on of A. The thin solid line is 
the exact result given by eq. (1.10) for the NN XX Z mo del. 



The broken line indicates its asymptote given by eq. (4.12) for 
large A. The open triangles (circles) show the data of Eq (Eq) 
extrapolated from finite size data (N = 6-24) for the IS XXZ 
model. The thick lines are guides to the eye. 



Fig. 5. Size dependence of the energy gaps (a) Eq and (b) Eq 
for the IS XXZ model with various values of A. The data with 
each A are fitted by the second-order polynomial. 



Fig. 6. S zz (q,uj) [= S xx (q,oj)] for the IS exchange model with 
A = 1 and N = 24. The intensity of each pole is proportional 
to the area of the circle. 



Fig. 7. The oj-dependence of S xx (tt, uj) for the IS exchange model 
with N = 16-24 and A = 1,2 and 5. The circles, triangles, 
squares, inverted triangles and diamonds are the data for N = 
16, 18, 20, 22 and 24, respectively. These data are obtained by 
processing of the poles with main intensity (all poles only for 
A = 1) via the LSD method. The inset shows the aj-dependence 
of NM xx (tt,uj) for A = 2. 
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Fig. 8. S xx (q, ui) of the IS exchange model with N = 24 and (a) 
A = 2, (b) A = 5 and (c) A = 10. The intensity of each pole is 
proportional to the area of the circle. 



Fig. 9. The oj-dcpcndcncc of S xx (tt,u>) for the NN exchange 
model with N = 16-24 and A = 1, 2 and 5. The symbols are the 
same as in Fig. 7. These data are obtained by processing of the 
poles with main intensity via the LSD method. The inset shows 
the w-dependence of NM XX (ir, u) for A = 2. 



Fig. 16. Excitation spectrum in the IS exchange model for N = 
100, e = 0.1 and SfL* = 1. The closed diam onds show the 
solutions obtained from the eigenvalue equation (5.13) using the 
two-soliton approximation. The solid line s are the lower and 
upper edges of the continuum given by eq. (5.17). 



Fig. 10. Dynamical quantities for the IS exchange model with 
A = 10. The symbols are the same as in Fig. 7. (a) The in- 
dependence of M xx (n,u>) and NM xx (it,lli) (inset) for N = 16- 
24. (b) The oj-dependence of S xx (n, ui). The data are ob- 
tained by processing of the poles except highest one via the LSD 
method. 



Fig. 11. The (^-dependence of M xx (tt,u>) and NM xx (tt, u>) (in- 
set) for the IS exchange model with N = 16-24 and (a) A = 6 
and (b) A = 7. The symbols are the same as in Fig. 7. 



Fig. 12. S zz (q, u>) of the IS exchange model with N = 24 and (a) 
A = 2 and (b) A = 5. The intensity of each pole is proportional 
to the area of the circle. 



Fig. 13. S xx (q) of the IS exchange model with N = 16-24 and 
A = 1, 2 and 5. The circles, triangles, squares, inverted triangles 
and diamonds are the data for N = 16, 18, 20, 22 and 24, respec- 
tively. The solid line is th e exact result for A = 1 (i.e., the HS 
model) given by eq. ( |4.1e| ). The inset shows the magnification 
in the region of small q. 



Fig. 14. S zz {q) of the IS exchange model with N = 16-24 and 
A = 1, 2 and 5. The symbols are the same as in Fig. 13. The 
solid li ne is the exact result for A = 1 (i.e., the HS model) given 
by eq. (4.19). The inset shows the magnification in the region of 
small q. 



Fig. 15. Excitation spectra in the IS exchange model for N = 24, 
e = 0.1 (i.e., A = 10) and S z ot = 1. (a) The closed dia monds 
show the solutions obtained from the eigenvalue equation (5.13) 
using the two-soliton approximation. The solid lin es ar c the 
lower and upper edges of the continuum given by eq. (5.17). (b) 
The pole positions by the exact diagonalization and the recursion 
method [see Fig. 8(c)]. The closed diamonds show the raw data 
in the region [0, w], and the open diamonds exhibit the data 
shifted from [tt, 2tt] to [0, it]. The plus indicates the pole position 
with very weak intensity (~ 10 — 6 ) compared to other ones. 



